Formation Scenario for Wide and Close Binary Systems 



Masahiro N. Machida^, Kohji Tomisaka^, Tomoaki Matsumoto^, and Shu-ichiro Inutsuka^ 



ABSTRACT 

Fragmentation and binary formation processes are studied using three- 
dimensional resistive MHD nested grid simulations. Starting with a Bonnor- 
Ebert isothermal cloud rotating in a uniform magnetic field, we calculate the 
cloud evolution from the molecular cloud core (n = 10'^ cm~^) to the stellar 
core (n ~ lO^^cm"'^), where n denotes the central density. We calculated 147 
models with different initial magnetic, rotational, and thermal energies, and 
the amplitudes of the non-axisymmetric perturbation. In a collapsing cloud, 
fragmentation is mainly controlled by the initial ratio of the rotational to the 
magnetic energy, regardless of the initial thermal energy and amplitude of the 
non-axisymmetric perturbation. When the clouds have large rotational ener- 
gies in relation to magnetic energies, fragmentation occurs in the low-density 
evolution phase (lO^^cm"^ ^ ^ ^ 10^^ cm~^) with separations of 3-300 AU. 
Fragments that appeared in this phase are expected to evolve into wide binary 
systems. On the other hand, fragmentation does not occur in the low-density 
evolution phase, when initial clouds have large magnetic energies in relation 
to the rotational energies. In these clouds, fragmentation only occurs in the 
high-density evolution phase {n > 10^''cm~^) after the clouds experience signifi- 
cant reduction of the magnetic field owing to Ohmic dissipation in the period of 
10^^ cm~^ < 10^^ cm~^. Fragments appearing in this phase have separations 
of < 0.3 AU, and are expected to evolve into close binary systems. As a result, we 
found two typical fragmentation epochs, which cause different stellar separations. 
Although these typical separations are disturbed in the subsequent gas accretion 
phase, we might be able to observe two peaks of binary separations in extremely 
young stellar groups. 
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Subject headings: binaries: general — ISM: clouds: ISM: magnetic fields — MHD — 
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INTRODUCTION 



Stars are born as binary or multiple systems. Observations have shown that, although 
the multiplicity depends on the stellar mass and surrou nding environment, about 60-80 % of 



field s tars are members of binary or multiple systerns fjHeintzl Il969l : lAbtfc Levyl Il976l : lAbt 
1983 : Duquennoy &: Mayor 1991; Fischer fc Marcy 19921 ). Thus, it has been established 



that a majority of main-sequence stars consist of binary or multiple systems. Between mul- 
tiples, the fraction o f binary (b), triple (t), and higher-order multiple (h) syste ms is b: t: h 



0.75: 0.175: 0.04 ( jPuquennoy &: Mayor Ill99ll : iTokovinin &: Smekhov 



cates that binary systems account for a large pe rcentage of the se systems ( iGoodwin et al. 



200211 ■ which indi- 



20071 ). Observations of star- forming regions (e.g.. lMathieulll994j ) have shown that the multi- 
plicity of pre-main sequence stars is larger than that of main-sequence stars. For example, in 
the Taurus star-forming region, it is expecte d that almost 100% of pre-main sequence star s 



are members of binary or multiple systems (iLeinert et al. 1 119931 : iKohler fc Leinert Ill998l ). 
Recently, extremely young protostars (i.e., Class prot ostars) have been observable with 



radio interferometer ( Loonev et al. 200ol: Reipurth 200ol ) and wide- field near- infrared cam- 



eras (IHaisch et al. 11200 



ii 



2004 



Duchene et al. 



20041 ). These observations have sliown_tliat 



stars already have a high multiplicity at the moment of their birth. In contrast, iLadal (120061 ) 
pointed out that the multiplicity of low mass main-sequence stars (< 0.5 Mq) is less than 
~ 50%. The rnultipli city of the field stars tends to decrease with decreasing stellar mass. 



Duchene et al. I (120061 ) denoted that the multiple system decays by dynamical disruptive in- 
teraction in a short timescale (< IMyr), and these systems evolve into an ejected single star 
or stable binary system. This mechanism is effective in multiple systems composed of lower 
mass stars, because the binding energy of these systems is small. Thus, it is considered that 
stars are born as binary or multiple systems, and some systems disrupt into single stars as 
time goes on. For this reason, binary and multiple frequencies decrease with time, especially 
in low-mass systems. 

Distribution of binary separation is known to be very wide (0.01 AU< r < lO ^AU) and 
flat, usually modeled as a lognormal with a mean ~ 30 AU (IGoodwin et al. 1120071 ). For field 
main-sequence stars, roughly a third of all companions have close separation (r < lAU). In 
addition, some binary systems have rotation periods of P < 10 days, which correspon ds to 
r < 0.1 AU of the separation when the total mass of 1 Mq is assumed (jMathieulll994j ). and 
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they are classif ied as close binary 0. These are the same for stars in the star-forming region. 
Mathieul (Il994j ) showed that the binary separation extends from 0.01 AU up to 1000 AU also 
in pre-main sequence stars. This indicates that both wide and close binary systems are 
formed in the star-formation process. However, the mechanism that determines the binary 
frequency and separation is still unknown. The frequency and separation distributions of 
binary are important clues to understanding the formation process of such stars. It is, 
however, difficult to observe the formation process of binary systems directly, because these 
systems are embedded in dense gas clouds, and star (or binary) formation timescales are 
expected to be extremely short. Thus, numerical calculations are necessary to understand 
the binary formation process including binary frequency and separation distributions. 



Under the spherical symmetry, iLarsonl (119691 ) , iTohlind (119821 ) , and lMasunaga fc Inutsuka 



(I2OOOI ) numerically investigated the evolution of a cloud from the molecular cloud core 
stage to that of protostar formation including the detailed radiative transfer. They found 
that the gas cloud collapses self-similarly to reach stellar density. According to the ther- 
mal evolution derived in their studies, cloud evolution is divided into four phases: the 



isothermal 



(10^ 



cm 



< 



< 10 



11 



cm 



adiabatic (10^^ cm ^ ^ n, 
, and protostar formation phases 



< 10^1 cm- 



c ^ 10^^ cm ^), second collapse 



> lO^^cm-^l 



where 



ric means the number density at the center of the cloud. The cloud collapses isothermally 
for Uc ^ 10^^ cm~^ (the isothermal phase or first collapse phase), then the gas around the 
center of the cloud becomes optically thick and collapses adiabatically after the central den- 
sity exceeds ric > 10^^ cm^'^ (the adiabatic phase). In the adiabatic phase, the adiabatic 
core (or the first core) surrounded by a shock is formed. When the central density reaches 



rir 



10 cm"-^, the molecular hydrogen begins to dissociate and the cloud collapses rapidly 
(the second collapse phase). After the molecular hydrogen is completely dissociated, the gas 
collapses adiabatically again and the protostar (or second core) is formed at ric — 10^^ cm~^. 
In these spherically symmetric calculations, however, cloud fragmentation or binary for- 
mation cannot be investigated, because these calculations only include the isotropic forces 
of the gravity and thermal pressure. Cloud fragmentation is controlled by the anisotropic 
forces caused by the cloud rotation, magnetic field, and turbulence. Thus, to investigate the 
fragmentation of the cloud, three-dimensional calculations are needed. 



The evol ution of rotating clouds in three dimensions hav e been investigat e d by 



authors (c.f.. iBodenheimer et al. II2OOOI : iGoodwin et al. 1120071 ). iMiyama et al. I (119841 ). and 



many 



Tsuribe fc Inutsukal (119991 ) calculated the evolutions of spherical clouds with initially uniform 
density and rigid-body rotation law in the isothermal regime. They have shown that such 



^ In this paper, for convenience, we define the binary system with separation of r < 0.1 AU as close 
binary, while the binary system with separation of r > 0.1 AU as wide binary. 



-4- 



clouds fragment if ao < 0.12 — 0.15, where ao = E'j^/E^ is the initial ratio of the thermal 
to the gravitational energy, and (3q = E-^/E^ is the initial ratio of the rotational to the 
gravitational energy. When fragmentation occurs in the isothermal phase, the separation of 
the fragments has a scale of 100 — 10^ AU, which corresponds to the Jeans length in this phase. 
However, observations have shown that molecular clouds are centrally co ndensed, and often 



exhibit density profil es resembling those of the Bonnor-Ebert sphere (e.g.. lAlves et al. 112001 



Kandori et al. 



20051 ). IBossI (119931 ) calculated the evolution of such clouds with exponential 



density profiles in the isothermal re gime, and he fo u nd th at fragmentation occurs only in 



highly unstable clouds of < 0.3. iKandori et al. I (120051 ) found, from t heir near-infrared 



i magin g observations, that many starless cores are near the stable state. iTachihara et al 



(I2OO2I ) observed 179 cores in the star-forming regions, and showed that many of the cores are 
nearly in virial equilibrium. Namely, their observations indicate that is not so small, but 
is close to ao — 1 B For this reason, it is expected that fragmentation rarely occurs in the 
isothermal phase. The evolution of rotating clouds with centr ally condensed density profiles 
fror n the isothermal to t h e adia batic phase are investigated by lMatsumoto fc Hanawal ( 20031 ) 
and ICha fc Whitworth I (120031). Thev found that clo ud rotation promotes fragmentation in 
the adiabatic phase. iMatsumoto fc Hanawal (120031 ) showed that fragmentation occurs in 
the adiabatic phase when the cloud core ha s a rotation energy oi 00 > 3 x 10 ~^, which 



i s suffi ciently smaller than the observations (iGoodman et al. I Il993l : ICasellil |2002| ). ICaselli 



(I2OO2I ) showed that the molecular clouds have 10~^ ^ /^o ^ 0.02 with the mean value of 
/So ~ 0.02. Thus, their results indicate that the major fraction of molecular clouds might 
fragment in the adiabatic phase. 

Many observations have shown that cloud cores have a significant non-thermal motion, 
which can be attributed to turbulence. The turbulence is thought to dissipate as a cloud 
collapses (e.g., lLarsonlll98l[ ). In the very early phase of cloud collapse, however, the cloud 
acquires the angular momentum from turbulence, besides the large-scale ordered rotation 
motion of the host cloud. Thus, turbulence promotes fragmentation, as for the cloud rotation. 
Numerical studies have sh own that more fragme nts appear in a cloud with larger turbulent 



energy in the initial state (iGoodwin et al. 1120071 ). The cloud evolution with larger turbulent 



energy, in wh ich an initial cloud has tu r bulent energy co r npara b le to the gravitational, is in- 
vestigated by ISate et al. I J2002I . l2003h , ISate fc Bound] J2005h , IPekado-Donate I J2004al lb[) . 
In their calculations, fragmentation frequently occurs in an isolated cloud core, and many 
binary or multiple systems appear. iGoodwin et "al~l J2004al lbl) calculated the evolution of 
clouds with smaller turbulent energy. They found fragmentation occurs and a binary system 
is born when the initial cloud has ~ 5% turbulence energy exceeding gravitational energy. 



^ The critical Bonnor Ebcrt sphere has ao — 0.7 ( Matsumoto fc Hanawa 2003h . 
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They have suggested that almost all stars are born as binary or multiple systems, even if a 
cloud core has a low level of turbulence. However, binary frequencies found in these turbu- 
lence calculations are higher than in observations. Since the magnetic fields were ignored in 
these calculations, these high binary frequencies may be caused by neglect of the magnetic 
effect. 

While cloud rotation and turbulence promote fragmentation, the magnetic field sup- 
presses fragmentation. The e volution of magnetized clouds f rom the isotherrnal to the adia - 
batic phase is i nvestigated bv iHosking fc Whitworth j2004h . |Machida et al. I J2004j . boOSbh . 
Ziegler ( 2005 ). Fromang et al. I ( 2006 ). and Price fc Bate! ( 2007 ). In their studies, fragmen- 
tation rarely occurs in strongly magnetized clouds. This is because the angular momentum, 
which promotes fr agmentation, is r emoved by the magnetic effect (e.g., magnetic braking, 
outflow, and jets). iMachida et al. I (l2005bl ) showed that fragmentation depends only on the 
ratio of angular rotation speed to the magnetic field strength of the initial clouds, and hence 
it is necessary for strongly magnetized clouds to rotate rapidly for fragmentation. 

Although fragmentation processes in the isothermal and adiabatic phase (n < 10^^ cm^^) 
are investigated by many authors, there are few studies of fragmentation in the second col- 
lapse and protostellar phases {n > 10^^ cm~^). When fragmentation occurs in the isothermal 
or adiabatic phases {n < lO^^cm"^), fragments have separations of 10 — 10^ AU, which cor- 
respond to the Jeans length in these phases. Fragments formed in these phases are expected 
to evolve into wide binary systems, when the binary systems maintain separation at their 
formation epoc h. However, about 20-30% of observed binaries have separations of < 10 AU 
( lMathieulll994l ). In the collapsing cloud, the Jeans length becomes < 10 AU after the cen- 
tral density exceeds n > lO^^cm"^. Thus, fragments with the separation of < 10 A U and 



Bate 



thus c lose b inary systems are ex p ected when fragmentation occurs for n > 10^^ cm~^ 
(119981 ) and IWhitehouse fc Batd (120061 ) investigated the evolution of unmagnetized clouds 
from 10^ cm~^ to 10^^ cm~^. According to their calculations, fragmentation does not occur 
in the high-density region, because the non-axisymmetric structure which appears in the 
adiabatic stage effectively removes angular momentum from the center of the cloud. Thus, 
they concluded that fragmentation does not occu r for n > lO^^cni " ^. Ho wever, they inves- 



tigated the cloud evolution in a few parameters. iGoodwin et al. I (120071 ) pointed out that 
cloud evolution in high-density gas is complicated, hence a large number of calculations and 
a statistical approach are needed to understand the fragmentation process. 



Banerjee fc Pudritzl (120061 ) investigated the cloud evolution of magnetized clouds for 



lO^cm"^ 10^°cm~^, and found that fragmentation occurs in the high density gas via 

a ring structure. They expected that this ring structure evolves into a close binary with the 
separation of Si?©. However, they adopted the ideal MHD approximation, which is not valid 
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in high-density gas as > 10^^ cm ^. A significant ma gnetic fiux is lost by O hmic dissipation 
in the density range of 10^^ cm~'^ ^ ^ 10^^ cm~'^ flNakano et ah 1 12002| ) and hence they 



overestimated the magnetic fiux in a collapsing cloud, especially in the high-density gas 
regions. In addition, they investigated the cloud evolution with only one parameter. 

In this study, we investigate the evolutions of magnetized cloud cores {uc — lO^cm"^, 
r = 4.6 X 10^ AU) until a protostar is formed {uc — lO^^cm"^, r ~ IRq) using three- 
dimensional resistive MHD nested grid simulations. We calculated 147 models with different 
magnetic field strengths, rotation speeds, and initial amplitudes of non-axisymmetric per- 
turbation. We found that the formation conditions for close and wide binary systems are 
related to the ratio of the rotational to the magnetic energy of the cloud core. 

The structure of the paper is as follows: The framework of our models is given in §2 
and in §3 the numerical method of our computations is shown. The numerical results are 
presented in §4. We discuss the fragmentation conditions in §5. 



MODEL 



Our initial settings are almost the same as those of Machida et al. ] J2006al lbl. [2007al lb[). 



To study the cloud evolution, we use a three-dimensional resistive MHD nested grid code. 
We solve the resistive MHD equations including the self-gravity: 



9p „ 
— + V 

dt 



p— + pyv ■ \/)v 



dt 



-VP 



dB 

'dt 



{pv) 
1 

V X (v X B)- 
VV = AnGp 



= 0, 

B X (V X B) - pV<^, 
]+riV^B, 



(1) 

(2) 

(3) 
(4) 



where p, v, P, B, r] and denote the density, velocity, pressure, magnetic fiux density, 
resistivity, and gravitational potential, respe ctively. The resistiv ity r] in equation Q is a 
function of the density and the temperature (INakano et al. II2002I ). We use the resistivity r] 



as the value a dopted inlMachida et a. 



calculated by iMasunaga fc Inutsuka 



, J2006bl . l2007j )7 To mimic the temperature evolution 
we adopt a piece-wise polytropic equation of 
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state as: 



P 



clo P 

Cs,0 Pc 
Cs,0 pc 



7/5 
7/5 
7/5 



1.1 



1.1 



Pc 

(rir 



5/3 



(P < Pc), 

(Pc < P < Pd), 

{pd< p< Pe), 

(P > Pe), 

, Pd = 



(5) 



where Cs,o = 190 ms~ , pc = 3.84 x 10~ "^g cm"'' {ric = 10 cm"''), pd = 3.84 x 10"^ g cm"'' 
{nd = 10^^ cm"^), and pe = 3.84 x 10"^g cm"'^ (rig = 10^^ cm"'^). For convenience, we 
define 'the protostar formation epoch' as that at which the central density {ric) reaches ric = 
10^^ cm"'^. We also call the period of Uc < 10^^ cm"^ 'the isothermal phase', the period of 
10^^ cm"'^ < He < 10^^ cm"'^ 'the adiabatic phase', the period of 10^^ cm"'^ < Uc < 10^^ cm"'^ 
'the second collapse phase', and the period of > 10^^ cm"'^ 'the protostellar phase.' 

As the initial condition of the cores, we use the density profile ii icreased by a factor 
/ (de nsity enhancement factor) from the critical Bonnor-Ebert sphere ( lEbertlll955l : iBonnor 
19561 ). Note that the critical Bonnor-Ebert sphere has a density contrast of 14 between the 
center and the surface of the cloud. The density enhancement factors are fixed as / = 1.68 in 
all models except for five models shown in §5.1. In §5.1, we adopt the density enhancement 
factors of /=1.4, 1.8, 2.8, 8.41, and 84.1. The density enhancement factor is related to the 
stability of the initial cloud. We discuss the relation of the density enhancement factor and 
the stability of the cloud in §5.1. For the critical Bonnor-Ebert sphere, the central density 
is equal to n^o = 10"^ cm"^ (po = 3.82 x 10"^°gcm"^), and the initial temperature is lOK. 
The radius of the critical Bonnor-Ebert sphere Rc = 6.45 Cs[47rG'pB£;(0)]"^''^ corresponds 
to Rc = 4.58 X 10^ AU. Outside this radius, we assume a uniform gas with a density of 
?t,be(-Rc) = 711 cm"^. The total mass contained in the critical Bonnor-Ebert sphere is 
Mc = 4.5 Mq, and our initial core is / times as massive as this value. 

Initially the cloud rotates rigidly Qq around the z-axis and has a uniform magnetic 
field Bq parallel to the 2;-axis (or rotation axis). The initial model is characterized by three 
nondimensional parameters: the magnetic-to-thermal pressure ratio. 



b 



47rpc,oC^,o ' 

the angular rotation velocity normalized by the free-fall timescale, 

^c,0 



UJ 



y/^nGp, 



(6) 



(7) 



'c,0 



and the initial amplitude of the non-axisymmetric perturbation A^p. We add m = 2- mode 
non-axisymmetric density perturbation to the spherical initial cores. Then, the density 
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profile of tlie core is denoted as 

p(r) = 



PBE(r) (1 + 6p) f for r < R, 
PBE(i?c)(l + M/ for r>R, 



where Pbe{i^) is tfie density distribution of the critical BE sphere, and Sp is the axisymmetric 
density perturbation. For the m = 2-mode, we chose 

Sp = A^{r/R,)cos2<P, (9) 

where represents the amplitude of the perturbation. 

We adopt & as & = — 4 as summarized in Table [H as = 0.003 — 0.2 as summarized 
in Table [21 and as A^p = 0.01, 0.2, and 0.4. The ratios of rotational (Pq) and magnetic (70) 
energies to the gravitational energy are summarized in each table. In addition, we estimate 
the mass-to-fiux ratio 

MM 

where M is the mass contained within the critical radius Rc and $ is the magnetic flux thread- 
ing the cloud. There exists a critical value of M/ $ below which a clo ud is supported against 



the gr avity by the magnetic field. For a cloud with uniform density, iMouschovias fc Spitzer 



(119761 ) derived a critical mass-to-fiux ratio 



where the constant C = 0.53 for uniform spheres {( = 0.48 by recent careful calculation of 
Tomisaka et"al~lll988al lb[). For convenience we use the mass-to-fiux ratio normalized by the 



critical value as 



/ norm \ / \ / cri 

The relation of parameter h and the mass-to-fiux ratio is also summarized in Table [H We 
made 147 models by combining the above parameters. 



NUMERICAL METHOD 



We adopt the nested grid method (for detail, see lMachida et al. Il2005al . l2006al ) to obtain 
high spatial resolution near the center. Each level of a rectangular grid has the same number 
of cells (= 128 X 128 x 8), although the cell width h{l) depends on the grid level /. The cell 
width is reduced by a factor of 1/2 as the grid level increases by 1 (/ / + 1). We assumed 
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mirror symmetry with respect to z=0. The highest level of grids changes dynamically. We 
begin our calculations with four grid levels (/ = 1,2,3, and 4). Box size of the initial 
finest grid / = 4 is chosen to be equal to 2Rc, where Rc denotes the radius of the critical 
Bonnor-Ebert sphere. The coarsest grid (/ = 1), therefore, has a box size equal to 2^ R^. A 
boundary condition is imposed at r = 2'* Rc , where the magnetic field and ambient gas rotate 
at an angular velocity of flo (for detail, see iMatsumoto fc Tomisakall2004j ). A new finer grid 
is generated whenever the minimum local Jeans length Aj becomes smaller than 8h{l^g,^), 
where h is the cell width. The maximum level of grids is restricted to /max = 30, in which 
the maximum density up to ric = 5 x 10^^ cm^^ can be calculated safely. Note that since 
the central density stopped increasing at ric < lO^^cm"^ in every model, the grids of / = 29, 
and 30 were not generated in our calculation. Since the densi ty is highest in t he fin est grid, 
the generation of a new grid ensures the Jeans condition of iTruelove et al. I (119971 ) with a 
margin of safe t y fact or of 2. We adopted the hyperbolic divergence B cleaning method of 



Dedner et al. I (120021 ). 



4. RESULTS 

In this section, we describe the evolution of the magnetized clouds parameterized by 
the magnetic field strength b, rotation rate uj, and initial amplitude of the non-axisymmetric 
perturbation A^. We compare the models with different b and u, but fixed in the following 
subsections. In §4.1, we show the evolutions for initial small non-axisymmetric perturbation 
A^p = 0.01. Then, the cloud evolutions for initial large non-axisymmetric perturbations are 
shown in §4.2 {A^ = 0.2) and 4.3 {A^ = 0.4). 



4.1. Model with Small Non-Axisymmetric Perturbations 

In this subsection, we show the evolution of clouds for a small non-axisymmetric per- 
turbation A^ = 0.01 at the initial state. Figure [1] shows the cloud structures around the 
center of the cloud for each model at the end of the calculation. In this figure, all the models 
have the same = 0.01, but different b and u. Each panel is placed according to the 
parameters b (x-axis) and u (y-axis). Models located in the upper- right region have strong 
magnetic fields and rapid rotations, while models located in the lower-left region have weak 
magnetic fields and slow rotations in the initial state. In the blue region, fragmentation 
occurs in the adiabatic (first core) phase (10^^ cm~^ ^ ""-c ^ lO^^cm"^), while fragmen- 
tation occurs after the second collapse [ric > 10^^ cm~'^) in the red region. On the other 
hand, fragmentation is never seen in the gray region. For convenience, we call the models 
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in which fragmentation occurs in the adiabatic phase (10^^ cm~^ < Uc < 10^^ cm~^) 'first 
fragmentation models', the models in which fragmentation occurs after the second collapse 
(tt-c > lO^^cm"^) 'second fragmentation models', and the models which never experience 
fragmentation 'non- fragmentation models.' 

In both the first and second fragmentation models, we stopped the calculations when 
the Jeans condition was violated outside the central deepest-level grid. In these models, 
after the fragments escape from the finest grid, the Jeans condition is violated in the coarser 
grid. In non-fragmentation models, we stopped the calculations in the protostellar phase 
[ric > 10^^ cm~^) after we confirmed that fragmentation was unlikely to occur. We had 
to stop the calculations in models near the lower-left [{b, u) = {0, 0.007), (0, 0.01)] and 
upper-right corners [{b, uj) = (4, 0.2)] (green region), before a protostar or fragmentation 
appeared. In the models with ti;) = (0, 0.007) and (0, 0.01), although we calculated the 
cloud evolutions for ~ 2000 yr after the first cores were formed, the formed first cores were 
stable and showed neither signature of collapse nor fragmentation. In addition, in these 
models, although we calculated the cloud evolutions for a long computational time of ~ 
1200 CPU-hours, these clouds do not collapse any more after the first core is formed. Note 
that it takes about 400 CPU-hours to calculate the cloud evolution until the protostar is 
formed for a typical model. The first cores in the models in the lower-left corner seem 



long-lived as shown in ISaigo &: Tomisakal (120061 ). On the other hand, the cloud indicates 



oscillation around the initial state without any indication of collapse in model {b, u) = {4, 
0.2). Since this model has the strongest magnetic field and most rapid rotation, the cloud 
does not collapse to form the first core. 

In the first fragmentation models (blue region), the cloud fragments into several pieces 
in the adiabatic phase (10^^ cm~^ ^ ""-c ^ 10^^ cm~^). Typical separations between two frag- 
ments are Rsep ~ 20 — 200 AU in these models. In general, these models have high rotational 
energies but low magnetic energies at the initial state. In the second fragmentation mod- 
els, the cloud fragments after the second collapse phase {ric > lO^^cm"^). The separations 
between the fragments in these models are Rsep ~ 0.01 — 0.5 AU. Compared with the first 
fragmentation models, the second fragmentation models have lower rotational energies and 
higher magnetic energies. In non-fragmentation models (grey region), protostars are formed 
without fragmentation. Compared with the fragmentation models, non-fragmentation mod- 
els have lower rotational energies and higher magnetic energies in the initial state. Figured] 
shows that fragmentation occurs in the lower density (or large scale; 20-200 AU) in the 
clouds having larger rotational energies and smaller magnetic energies, and fragmentation 
rarely occurs in clouds having smaller rotation energies or larger magnetic energies. In the 
models shown in Figure [H the clouds evolve keeping an axisymmetry because of the small 
non-axisymmetric perturbation. Thus, fragmentation is induced via ring structure in many 
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models, and the formed protostars have outer axisymmetric disks. 



4.1.1. Typical Model (A^ = 0.01) 



First Core Formation To see the evolution of a molecular cloud into a protostar in more 
detail, we plot a number of snapshots for a typical model in Figured The model corresponds 
to {b, uj, A^) = (0.1, 0.05, 0.01). As listed in Tables 1 and 2, this model has 70 = 9.59 x 10^^ 
and Po = 8.24 x 10~^, respectively. Figure [2]a shows the initial density distribution along the 
z = 0, y = 0, and x = planes. We adopted a Bonnor-Ebert density profile and the initial 
cloud has a spherical structure with the central density of ric = 1.68 x lO^cm"^. 

Figure[2f) shows the cloud structure when the central density reaches Uc = 7.3 x 10^^ cm~'^ 
(adiabatic stage). In this model, when the central density reaches Uc — 5.1 x lO^^cm"^, a 
shock appears (r ^ 2AU) and a first core is formed at the center of the cloud @ . The 
first core is represented by a thick red line in Figure The first core has an oblate shape 
as shown in the middle and lower panels of Figure [2]&. Since the cloud collapses along the 
z-axis, the central region becomes oblate before the first core formation. To evaluate the 
core shape, we define the oblateness as ^ob = {hihsY^"^ /h^, where hi, hs, and hz mean the 
major, minor, and 2;-axis length derived from the moment of inertia for the high-density gas 
of p > O.lpc according to Matsumoto & Hanawa (1999). Figure [3]a shows the evolution of 
this oblateness. In the isothermal phase, the oblateness continues to increase from the initial 
state and a disk-like structure is finally formed, because both the Lorentz and centrifugal 
forces make the cloud oblate. The oblateness reaches ^ob — 5, and it reaches a peak when 
the central density reaches ~ 10^^ cm~^. In the adiabatic phase, the oblateness begins to 
decrease because the central region collapses adiabatically. The oblateness is ^ob — 1.2 when 



rir ~ 



10^^ cm ^. 



After first core formation, outflow appears near the first core. In the middle and lower 
panels of Figure it is shown that the gas is outflowing from the central region inside 
the thick orange lines. After the gas is supported by the thermal pressure, the magnetic 
field lines are strongly twisted because the rotational timescale becomes shorter than the 
collapse timescale. The outflow is driven due to the twisted magneti c field lines and rotation 
of th e first core. Outflow s froni the first cores are showri also in Tomisaka _( 1998 . I2OO0I . 
2002yBangriee fc Pudrit3 J2006h . [Matsumoto fc Tomisakal J2004h . iFromang et al. I tood ). 



and 



Machida et al. 



20041 . l2005bl . l2007al jbl) 



In this paper, we use 'first core formation' when the core surrounded by a clear shock boundary is 
formed in the adiabatic phase. 
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Figure [2]c shows the snapshot around the first core at the age of 263.4 yr after the first 
core formation. Although the first core has an almost axisymmetric structure at its formation 
epoch (Fig. ^ upper panel), the non-axisymmetric perturbation grows and a spiral pattern 
appears inside the first core, as shown in the upper panel of Figure [2]c. To evaluate the 
degree of non-axisymmetric pattern, we define the axis ratio as £ar = hi/hg — 1, according 
to Matsumoto & Hanawa (1999). Figure [3f) shows the evolution of the axis ratio. The axis 
ratio decreases for 1 0^cm~^ < Ur 5 , lO^crn "^, because we start with a nearly equilibrium 
state (for detail, see iMachida et "aLl boOSal lbl) . Then, the axis ratio increases in proportion 
to £ar ~ P^/^ for 10^ cm ^ ^ "^c ^ 10^^ cm ^. This growth rate corresponds to the result 
of the linear theory. (The evolutioii of a xis ratio can be described by the linear theory 



developed by iHanawa fc Matsumoto Il2000l and lLaill2000l .) In the linear regime, the degree 
of the deformation (i.e., the amplitude of the axis ratio) grows as ^ar oc with n being a 
constant for a given polytropic index 7. For the isothermal gas (7 = 1), the power index 
of the density is n = 1/6 (^ar oc p^^^). The power index of the density n decreases with 
increasing polytropic index 7. For 7 > 1.1, the index n becomes negative. Namely, the axis 
ratio increases (n > 0) for 7 < 1.1, whi le the axis ratio decre ases (n < 0) for 7 > 1.1. The 
power index n is given in Figure 14 of lOmukai et al. I (120051 ). In the adiabatic phase, the 
axis ratio decreases for 10^^ cm~^ ^ "^c ^ lO^^cm"^. This is because the polytropic index 7 
is changed from 7 = 1 to 7 = 1.4. 

The axis ratio increases very rapidly near ric — 3 x 10^^ cm~^, because collapse of the gas 
slows at this phase. The non-axisy mmetric perturbation can grow sufficiently in the quasi- 



static d isk for the isothermal phase (INakamura fc Hanawalll997l : lMatsumoto &: Hanawalll999 



see also Machida et al. 2007c ). Durisen et al. ( 19861 ) shows that the bar structure develops 
with the bar mode instability, wher i the ratio o f the rotational to gravitational energies 
(/3) of the core exceeds P > 0.274. ISaigo et al. I (120071 ) found that the first core forms a 
bar or spiral pattern after exceeding /5 > 0.27. Figure [3]c shows the evolution of angular 
velocity normalized by the square root of the density [^^.q/ {AnGpY^'^ = tUc'- the normalized 
angular velocity]. This panel shows that the normalized angular velocity saturates around 
ujc — 0.1—0.2 in the isothermal phase, while it increases in the adiabatic phase. After the first 
core is formed, the normalized angular velocity oscillates intensely ar ound 0;^ 0.2 — . 6. The 
evolut i on of a ngular velocity in the collapsing cloud is described by iMachida et al. I (l2005al . 
2006al . l2007al ). Since /? is related to Uc hj P = for a spherical cloud with a constant 



density, the first core is expected to have /? = 0.04 — 0.36. Thus the bar mode i nstability can 
be induced in the first core by the dynamical instability (iDurisen et al. 1119861 ). As a result, 
the first core shown in Figure [2]c upper panel has a spiral structure, which is thought to be 
caused by nonmagnetic instability. 



Figure [2]c upper panel shows two density peaks near the center of cloud, which are 
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caused by the fragmentation of the first core. In this model, fragmentation and subsequent 
merger are repeated twice in the range of 10^^ cm~^ ^ ""^c ^ 10^^ cm"'^. The first fragmen- 
tation occurs at the density ric — 2 x 10^^ cm~^, and then, the fragments are merged in 
a short timescale. The second fragmentation occurs at nc ~ 8 x 10^^ cm~^, and then the 
density exceeds the H2 dissociation density and the second collapse occurs in each fragment. 
As shown in Figure [3]c, the normalized angular velocity decrease rapidly in the range of 
lO^^cm"'^ ^ '^c ^ lO^^cm"^. This is driven by the angular momentum redistribution into 
orbital and spin motions. The two fragments are merged to form a nearly spherical core 
at the density ric ~ 2 x 10^*^ cm~^. The normalized angular velocity increases again after 
reaching n,, ^ 2 x 10^*^ cm~'^, because the fragments bring the angular momentum into the 
center of the cloud again. After merger, the central region changes shape from spherical to 
bar- like. At the protostar formation epoch (nc = 10^^ cm~^), the central region has an axis 
ratio of £ar = 2.8 and a bar-like protostar is expected in this model. 

Protostar Formation Figures show the cloud evolution after the protostar is formed. 
The axis ratio continues to increase in the protostellar phase, and an elongated bar is formed 
around the center of the cloud. When the axis ratio reaches e^v = 10.2 at ric — 1.8 x 10^^ cm~^, 
the bar fragments into three pieces as shown in Figure [2]d. When the fragmentation occurs, 
the separation between fragments is equal to Rgep — 0.031 AU. After fragmentation, frag- 
ments approach each other, and they are merged into a single core as shown in Figure [2]e. 
After the merger, fragmentation occurs again through a ring-like structure as shown in Fig- 
ure [2]f. The final fragments survive separately, until we stop the calculation. At the end 
of the calculation, the separation between fragments is equal to Rsep — 0.03 AU, and each 
fragment has a mass of M = 3.1 x IO^^Mq. 



Magnetic Field Figure [3H shows the evolution of the magnetic field strength normalized 
by the square root of the central density \ B^,c./(4: 7 ic^ nPcY ^'^- the normalized magnetic field 
strength]. Since, in lMachida et al. I (l2005al . l2006al . l2007al ) . we already discussed in detail the 
evolution of the magnetic field in the collapsing cloud, in this paper we briefly show the evolu- 
tion of the magnetic field. In the isothermal phase, the normalized magnetic field strength is 



saturated at a constant level according to the magnetic fiux-spin relation of iMachida et al. 



(l2005al . l2006al ). then the magnetic fields are dissipated by Ohm ic d issipation in the range o f 



10^2 cm 



-3 < 



< 10^^ cm-^ as shown in lNakano et al. I fl2002h and lMachida et al. I fl2007al ). 



The rapid decrease in Figure [3]c? for 10^^ cm~^ ^ '^c ^ 10^'' cm~'^ corresponds to this Ohmic 
dissipation. The first core, which is traced by two spiral arms occupying r < 3AU, has a 
large beta /5 ~ 10 — 10^, because the gas density exceeds ric > 10^^ cm~^. Thus, in the 
adiabatic phase, the magnetic field barely affects the cloud evolution inside the first core. 
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On the other hand, the region outside the spiral arm contains a strong magnetic field of 
/3 ~ 1, because the gas is well coupled with the magnetic field for ric < lO^^cm"^. In the 
second collapse phase, since the Ohmic dissipation becomes ineffective for the gas in which 
the ionization degree is recovered by the electron for alkali metals, the magnetic field can be 
amplified again, as shown in Figure [3](i. At the end of the calculation, each protostar has a 
magnetic field of l.OSkG at its center. 

4.2. Model with Large Non-Axisymmetric Perturbations 

In this section, we show evolutions of clouds with large non-axisymmetric perturbations 
{A^p = 0.2 and 0.4) at the initial state. 

4.2.1. Models with = 0.2 

Figure m shows cloud structures around the center of the cloud at the end of calculations 
for models with A^ = 0.2. In the upper-left region (weak magnetic field and fast rotation), 
fragmentation occurs in the adiabatic phase, and two fragments have wide separations (20 < 
Rsep ^ 100 AU). On the other hand, fragmentation occurs in the second collapse or after the 
protostellar phase with narrow separations (-Rsep ^ 0.4 AU) in the lower region (red region). 
Clouds having larger magnetic energies and smaller rotational energies at the initial state 
lead to formation of a protostar without fragmentation (grey region). Thus, similarly to 
models with small non-axisymmetric perturbation A^ = 0.01, the magnetic field suppresses 
fragmentation, but the rotation promotes it also in the models with large non-axisjnumetric 
perturbation A^ = 0.2. However, three regions of first (blue region), second (red region), and 
non- fragmentation (grey region) are clearly divided in Figure [1], while the boundary between 
these three regions is complicated in Figure |H Figures [T] and H] show that, even though the 
initial clouds have the same magnetic field strength and rotation rate, fragmentation epoch 
depends on the initial amplitude of the non-axisymmetric perturbation A^. For models 
of (6, u) = (0, 0.07), (0, 0.05), (0.001, 0.05), and (0.01, 0.05), fragmentation occurs in 
the adiabatic phase in the models with A^p = 0.01 (Fig. [T]), while fragmentation occurs in 
the second collapse or protostar phases in the models with A^ = 0.2 (Fig. H]). The epoch 
of fragmentation is delayed with increasing perturbation amplitude A^. In addition, some 
models in which fragmentation is found in Figured] (models with A^ = 0.01) do not fragment 
when the initial cloud has a large non-axisymmetric perturbation A^p = 0.2. For models of 
(6, u) = (1, 0.03), (1, 0.05), (1, 0.07), and (1, 0.1), fragmentation occurs in the models with 
A^ = 0.01 (Fig. d]), while fragmentation does not occur in the models with A^ = 0.2 (Fig. H]). 
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As shown in iMachida et al. I (I2005bl . l2007bl ). there are two modes of fragmentation: ring 
and bar fragmentations. The fragmentation mode is related to the axis ratio and rotation rate 
of the central core. When the cloud has a small axis ratio and rapid rotation, fragmentation 
is induced via a ring-like structure. On the other hand, when the cloud has a large axis ratio 
(^ar ^ 10), fragmentation occurs via a bar-like structure, irrespective of the cloud rotation 
rate. Non-axisymmetric structure has two distinct effects on fragmentation: suppression and 
promotion of fragmentations. When a core (the first core or second core) has a moderate axis 
ratio (1 ^ £ar ^ 10), the angular momentum is removed from the core by gravitational torque 
due to the non-axisymmetric pattern and thus fragmentation is suppressed. In contrast, 
when the core has a large axis ratio (^ar ^ 10), fragmentation occurs through the bar 
instability, irrespective of the cloud rotation rate. Thus, compared with models having small 
non-axisymmetric perturbation = 0.01 (Fig. [T]), in models with large non-axisymmetric 
perturbation A^p = 0.2 (Fig. H]), fragmentation is suppressed in some models [e.g., models 
(6, u) = (0.001, 0.05), and (0.01, 0.05 )] by removal of angular momentum owing to the 
gravitational torque, although fragmentation is promoted in some models by the prominent 
bar mode [{b, u) = (0.001, 0.03)]. However, comparison between Figures [1] and H] indicates 
that only the models located near the boundaries between these three regions (first, second, 
and non-fragmentation) are influenced by the initial amplitude of the non-axisjnumetric 
perturbation. The trend that the magnetic field suppresses fragmentation while rotation 
promotes it is the same for both the models with A^ = 0.01 and A^ = 0.2. 



4.2.2. Typical Model (A^ = 0.2 ; 

Figure shows the cloud structure for the model {b, u, A^) = (0.01, 0.05, 0.2). In 
this figure, the structures around the center of the cloud with different scales are plotted 
for the final state of the calculation. The spatial scale of each successive panel differs by a 
factor of four and thus the scale between Figures [5]a and [5]f differs by a factor of 1024. In 
this model, fragmentation occurs in the second collapse phase, while fragmentation occurs 
in the adiabatic phase for the model with the same b and ut but smaller A^ = 0.01, as shown 
in Figure [TJ This is because the angular momentum in the model with A^p = 0.2 is more 
effectively transferred from the center of the cloud than in the model with = 0.01. This 
stabilizes the first core for fragmentation. 

Middle and lower panels of Figures Ela and b show a thin disk, which is formed in 
the adiabatic phase. Red-dotted lines in Figures [5]a and b indicate a shocked region which 
corresponds to the first core. Since this model has a large amount of non-axisymmetric 
perturbation A^p = 0.2 in the initial state, the non-axisymmetric pattern grows sufficiently 



- 16 - 



and a spiral appears just after the first core formation. After that, the central compact core 
is separated from the ambient ring-like structure. The upper panel of Figure [5]6 shows that 
the first core is composed of a central compact core and a surrounding ellipsoidal structure. 
There is a clear gap between the central core and ambient medium. In the middle and lower 
panels of Figures [5]a and b, orange lines indicate the outflow. Inside these lines, the gas is 
outflowing from the center of the cloud with ~ 1.2kms~^. The outflow appears when the 
central density reaches 2.1 x 10^^ cm~'^ in this model. 

Figures [5]c and d indicate that the central region has entered the second collapse phase 
because the density exceeds ric > lO^^cm"^, while the surrounding ring-like structure stays 
in the adiabatic phase. Although we cannot see any internal structure in the compact core in 
Figures Eh and b, Figures a nd d show the central region forms a spiral structure. Similar 
structure is seen in lBatd ( 1l998l ). who calculated the evolution of an unmagnetized cloud core 
having the parameter of u = 0.07 with an initial number density of n^o = 3.6 x 10^ cm~^. He 
found that the cloud collapses to form a single star without fragmentation after the second 
collapse phase, because the angular momentum is effectively removed from the center of the 
cloud by the spiral patterns. Although our model contains the magnetic field, the magnetic 
field becomes extremel y weak in the second co llapse phase due to the Ohmic dissipation. 



as shown in §4.1.1 and iMachida et al. I (l2007al ). thus, the magnetic field barely affects the 



cloud evolution in this phase. For this reason, both the unmagnetized and magnetized cloud 
are thought to follows a similar evolutional pattern in the second collapse phase. 

In the upper panels of Figures[5]e and /, four fragments appear inside the spiral structure. 
In this model, the first fragmentation occurs when the central density reaches ric = 1.1 x 
lO^^cm"'^, and two clumps appear. They are located at {x, y) ~ (±0.01 AU, ±0.07 AU) as 
shown in the upper panel of Figure [5]f. After that, fragmentation occurs again in each clump 
at ric — 4 x 10^^ cm~^, and as a result four fragments are formed as shown in Figure [5]f. 
Since the gas density of these four fragments exceeds > 10^^ cm^^, these fragments are 
protostars in our definition. Thus, we found a quadruple protostellar system. At the end 
of the calculation, the protostars have the furthermost separation of -Rsep ~ 0.15 AU. To 
investigate whether these protostars actually evolve into a quadruple stellar system without 
merger, further long calculation is needed. 

In the protostellar phase, jet appeared from the center of the cloud with the velocity 
of ~ 18kms~^. The boundary between the inflow and jet is represented by orange lines in 
the upper and middle panels of Figures [5](i and /, inside which the gas is outflowing from 
the protostar. In this model, two flows appear: a low velocity flow is driven by the rotation 
of the first core, and a high velocity flow is driven by the spinning motion of the protostar. 
Since we investigated in detail the evolution and mechanism of outflow and jet in a separate 
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paper (IMachida et al. Il2007bl ). we do not describe these further here. The outflow and jet 



driven from the binary system will be investigated in a subsequent paper. 



4.2.3. Models with = 0.4 

Figure E] shows the cloud structures around the center of the cloud at the end of the 
calculation for models with A^p = 0.4. The fragmentation properties (fragmentation epoch, 
scale, number of fragments, and structure) shown in this figure are almost the same as those 
in Figures [1] and HI In the following, we only describe the differences between the models 
= 0.4 (Fig. ED, = 0.01 (Fig. OP) and A^ = 0.2 (Figgl). 

Firstly, we compare the models having the same parameters of {b, u) = (0.001, 0.07), 
which are located in the third row and the second column in Figures [H HI and [61 Fragmenta- 
tion occurs in the adiabatic phase with models having A^p = 0.01 and 0.2, while no fragmen- 
tation occurs through any of the phases of cloud evolution in the model with A^ = 0.4. In 
Figure O for this model, we can see that the central region around the protos t ar ind icates a 



bar-like structure in a large scale, and a spiral structure in a small scale. iBatd (Il998l ) showed 
that the spiral structure effectively removes the angular momentum from the center of the 
cloud, and a single protostar is directly formed without fragmentation. Also in this model, 
since the angular momentum is effectively removed owing to the bar or spiral structure, the 
protostar is thought to form without fragmentation. 

Next, we describe the models with (6, u) = (0.001, 0.05), and (0.01, 0.05), which are 
located at the fourth row and the second and third columns. Figures [H IH and M show that 
fragmentation occurs in the adiabatic phase with wide separations (~ 40 AU) for the models 
with A^p = 0.01 and 0.4, while fragmentation occurs with a narrow separation (~ 0.1 AU) 
in the second collapse or protostellar phase for the models with A^ = 0.2. According 
to the degree of non-axisymmetric pattern (i.e., the axis ratio), these models experience 
three different types of evolutions: (i) if the core has a small axis ratio, fragmentation 
occurs via a ring-like structure in the adiabatic phase owing to the rapid rotation {A^ = 
0.01), (ii) if the core has a moderate axis ratio, fragmentation occurs in the second collapse 
phase, because the angular momentum is effectively removed by the bar-like structure in 
the adiabatic phase {A^p = 0.2), and (iii) if the core has a large axis ratio, fragmentation 
occurs via a bar-like structure in the adiabatic phase, because a considerably elongated bar 
is formed in the adiabatic phase and such a structure is unstable {A^ = 0.4). Therefore, the 
separation between fragments does not represent the initial amplitude of the nonaxisymetric 
perturbations. 
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Thirdly, we compare the models having the same [b, uj) = (0.001, 0.03), which are 
located in the fifth row and the second column in Figures [T|, HI and [61 In these models, 
fragmentation occurs in the second collapse phase for the model with = 0.01, while 
fragmentation occurs in the adiabatic phase for the models with = 0.2 and 0.4. This 
difference is understood as follows: For the model with a small non-axisymmetric perturba- 
tion {A^ = 0.01), since both the rotation rate and the non-axisymmetric perturbation are 
small when the first core is formed, neither ring nor bar fragmentation occurs. On the other 
hand, fragmentation occurs in the bar for models with larger non-axisymmetric perturba- 
tions {A^p = 0.2 and 0.4), because a considerably elongated bar is formed at the first core 
formation epoch. 

Finally, the strong magnetized models b > 1 are described. Figures [U HI and M clearly 
show that fragmentation becomes harder with increasing A^^ especially in the strong mag- 
netized models. For models with A^ = 0.01 (Fig. [1]), fragmentation occurs in the second 
collapse phase for models having u = 0.03, 0.05, 0.07 and 0.1 for b = 1, and 07=0.07 and 0.1 
for b=4. In contrast, for models with A^ = 0.4 (Fig. H]), no fragmentation occurs in models 
having the same b and u. The increase of th e angular velocity is sup pressed by the magnetic 



effect (e.g., magnetic braking and outflow; iMachida et al. I l2007al ) in strongly magnetized 



models, and then these models manage to fragment in the second collapse or protostellar 
phase when the initial cloud has a small amount of the non-axisymmetric perturbation of 
A^p = 0.01. In contrast, for a models with large A^, since the angular momentum is removed 
not only by the magnetic effect but also the gravitational torque as shown in Figure Elc-/, 
fragmentation is more severely suppressed. 

Figures [1], HJ and [6] show that the magnetic fleld strength, angular velocity, and non- 
axisymmetric perturbation play an important role in fragmentation, and they are closely 
related. Thus, it is rather hard to derive the fragmentation condition from small numbers 
of calculations. 



4.2.4. Typical Model (A^ = 0.4) 

Figure [7] shows the cloud structure on the z = plane for the model of {b, u, A^) 
= (0.01, 0.07, 0.4). In this flgure, the structure around the center of the cloud is plotted 
with different scales at the end of the calculation. The spatial scale of each successive panel 
differs by a factor of eight and thus the scale between Figures [7]a and d differs by a factor 
of 512. In this model, the flrst core is formed at nc ^ 3 x lO^^cm"^. The flrst core repeats 
fragmentation and merger in the adiabatic phase. As shown in Figure [7]a, in a large scale, 
the high-density region is composed of a central compact core and a surrounding ring-like 
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structure. At the end of the calculation, the central compact core has the maximum density 
of ric = 1.1 X lO^^cm"^, and a mass of M = 0.065 M©. The second collapse occurs in the 
central core at 459 yr after the first core-formation epoch. Then, the central region fragments 
again at n,, = 8 x lO^^cm"^, and two protostars are formed, as shown in Figures [7]c and 
d. Two fragments (or protostars) in Figure [7]d have a mass of M = 4.4 x 10~^ Mq, and a 
separation of i?sep = 0.073 AU at the end of the calculation. 

On the other hand, in the surrounding ring-like structure, there are two clumps at (x, 
y) ~ (±5 AU, =f25 AU) as shown in Figure[7]a. These clumps have masses of M ~ 0.036 M©. 
The clumps have the maximum density of tt-c = 8 x 10^^ cm~^ and they continue to collapse 
until the calculation ends. Thus, these clumps are expected to be more compact as time 
goes on, and then, a single or several protostars are formed in each clump. 

As shown in Figure [3, four clumps appear in this model. Clumps formed in the adiabatic 
phase have a wide separation of i?sep — 26 AU, and they seem to be a triple system in a 
large scale (Fig. [7]a). In the central compact core, there are two fragments, each of which 
contains a close binary system. This system is thought to evolve into a hierarchical star 
system, which is composed of two protostars with a large separation and two binary systems 
with a narrow separation, if fragments do not merge with each other and separations remain 
almost unchanged. 



5. DISCUSSION 

5.1. Initial Cloud Stability and Its Evolution 

In the gravitationally contracting core, the thermal pressure, rotation, and magnetic 
field work against the self-gravity. Generally, in the initial clouds, the former three forces are 
characterized by the ratios to the gravitational energy, (thermal energy), (3q (rotational 
energy), and 70 (magnetic energy), and the clouds are specified by the combination of these 
three parameters. In the previous section, we investigated the evolutions of clouds with the 
same ao (ao = 0.5) but different Pq and 70. The magnetic field and rotation are important 
for cloud fragmentation, because these forces cause the anisotropic patterns which promote 
fragmentation. However, the initial ratio of the thermal to the gravitational energy, which is 
related to the stability of the initial cloud, may affect cloud fragmentation. In this section, 
we compare the evolutions of clouds with the same f3o and 70 but different oq. 

We investigate the evolutions of clouds with five different (oq =0.01, 0.1, 0.3, 0.5, 
and 0.7). The cloud with ao = 0.01, which has 1% of the thermal energy to the gravitational 
energy at the initial state, is highly gravitationally unstable, while that with ao = 0.7 is 
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initially nearly in equilibrium. We fixed the three other parameters as [b, uj, A^^,) = (0.01, 
0.05, 0.01). Figure [8] shows the evolution of the oblateness ^ob (Fig. [8]a), axis ratio ^ar 
(Fig. [S]?'), normalized magnetic field strength B^c/ {'^'^c^oPcY^'^ (Fig- Elc), and normalized 
angular velocity Q^J{4nGp^y/^ (Fig. EH). 

Firstly, we show the evolution of clouds with large ao (ao = 0.3, 0.5 and 0.7). These 
models are initially more stable. In Figure [H these models show similar evolutions for 
oblateness, normalized magnetic field strength, and normalized angular velocity. After the 
central densities exceed Uc ^ 10^^ cm~'^, the central regions have almost the same degrees of 
oblateness, magnetic fields, and angular velocities as shown in Figures [8]a, c, and d. Shock 
fronts appear and the first cores surrounded by the shock waves are formed at ric — 10^^ cm~^ 
in these models. The first cores have almost the same properties except for the axis ratios. 
Comparing all five models, the formed first cores are sho wn t o have lar ger axis ratios in 
the models with smaller ao- iHanawa fc Matsumoto I (120001 ) and iLail (120001 ) showed that the 
growth rate of the non-axisymmetric pertu rbation is larger in a more unstable cloud core 
from their linear analysis (see also Fig. 14 of lOmukai et al. II2005I ). Since the cloud evolution 
converges to a se lf-similar solution as the cloud collapses even when the cloud is rotating 
and magnetized (IMachida et al. I l2005al . l2006al ) , the growth rate of the non-axisymmetric 
perturbation is also converged in the collapsing cloud. Thus, the difference in the axis ratio 
which reflects the amplitudes of the non-axisymmetric perturbation comes from the growth 
in the early phase of the contraction. For this reason, the initially more unstable cloud (i.e., 
the cloud with small ao) has a large amplitude of non-axisymmetric perturbation. 

Figure [9] shows the cloud structures near the center of the cloud just after fragmenta- 
tion. Fragmentation occurs via a ring-like structure in the models ao = 0.5 and 0.7, while 
fragmentation occurs without the ring formation in the model with ao = 0.3. This difference 
in the fragmentation pattern comes from the difference in the axis ratio of each model: the 
model with ao = 0.3 has a larger amplitude of non-axisymmetric perturbation than those of 
models with ao = 0.5 and 0.7 at the first core formation epoch. However, the scale and frag- 
mentation epoch are almost the same for all models with ao > 0.3, as shown in Figures [9]6-(i. 
Thus, even though the fragmentation patterns are different, the fragmentation scale and its 
epoch are hardly changed for models ao > 0.3. 

Next is shown the evolution of the highly unstable cores with ao = 0.1 and 0.01. Figure[8] 
indicates that the evolutions of oblateness, angular velocity and magnetic field of the models 
with ao = 0.1 and 0.01 differ appreciably from those of models with ao > 0.3. In addition, 
the first core in the model of ao = 0.1 is considerably smaller than those in models with 
larger ao (Fig. [9]). In the model with ao = 0.01, an extremely thin d isk with ^ob — 50 is 
formed in the isothermal phase (Fig. [8]a). iTsuribe fc Inutsukal (119991 ) expected these thin 



- 21 - 



disks to fragment in the isothermal phase when the cloud collapses isothermally for a long 
time. In our calculation, however, fragmentation does not occur in the isothermal phase, 
because the central region begins to behave adiabatically just after the oblateness reaches 
its peak. Then, an elongated bar is formed in the adiabatic phase, and the bar fragments 
at the central density Uc — 10^^ cm~^. The fragmentation scale of the model of ao = 0.01 is 
~ 2AU, which is considerably smaller than in models with > 0.3. Thus, fragmentation 
size becomes smaller with decreasing ao for < 0.3. 

In summary, for models with ao > 0.3, the clouds show similar evolutions, and the first 
cores have almost the same properties. On the other hand, clouds with ao < 0.3 evolve in a 
completely different way, and the first cores have different properties. However, observations 
have shown t hat molecular cloud co res are observed in nearly thermal equilibrium against 



gravity (e.g., iTachihara et al. II2002I ). which indicates ao ~ 1. As a result, the evolution of 
molecular cloud cores is specified by only (3o and 70 since clouds with ao > 0.3 have the same 
evolutional properties when the clouds have the same /?o and 70. 



5.2. Scales and Epochs of Fragmentation 



We calculated 147 models with different sets of initial cloud rotation, magnetic field, 
thermal energy, and amplitude of the non-axisymmetric perturbation. In 102 out of 147 
models, we observed fragmentation. To investigate the scales and epochs of fragmentation 
statistically, for all models that show fragmentation in Figures [1], HJ and [6l the number 
densities at the fragmentation epoch, and the farthermost separations between fragments at 
the end of the calculation are plotted in Figure dU] upper left panel (main panel). Right and 
lower panels in Figure [TD] show the histogram of fragmentation scales and epochs measured 
in the central density, respectively. Although we calculated the cloud evolution from ric = 
10^ cm~^, we observe no fragmentation in the isothermal phase (ric < 10^^ cm~^). Thus, only 
the range of 11^ > 10^ cm~^ is plotted. For convenience, we divide fragmentation models into 
three groups (Groups A, B, and C) by the evolution stage when fragmentation occurs. Models 
in Groups A, B, and C, respectively, show fragmentation in the adiabatic (10^^ cm^'^ < Uc < 
10^^ cm~'^; A), second collapse (10^^ cm~^ < Uc < 10^^ cm~^; B), and protostellar phase {uc > 
10^^ cm~'^; C). The dotted line in the main panel shows the Jeans length, which is derived 
from the relation of the temperature and density assumed in our calculation. Figure [TUl 
shows that points are distributed near the dotted line, which indicates that fragments have 
a separation nearly equal to the Jeans scale when they are born. Note that models are 
distributed slightly above the Jeans length in Figure [T0|, because a larger critical wavelength 
is expected for rotating magnetized gas, while the Jeans length plotted here is expected to 
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be a spherical symmetry. 

Figure [TU] shows that there are two distinct fragmentation epochs, namely ~ 10^^ cm~^ 
and ric ^ lO^^cm"^. The lower panel shows models belonging to groups B and C to be 
smoothly distributed, but there is a clear gap between groups A and B for 10^^ cm^^ < 
ric < lO^^cm"'^. Models in group A are distributed only in the range of 10^^ cm~^ < 
Tic < 10^^ cm~'^, although the adiabatic phase lasts for 10^^ cm~'^ < ric < 10^^ cm~^. In 
group A, the fragmentation epochs are bunched into the first core formation epoch. The 
first core is formed at lower density when the host cloud has a strong magnetic field or 
rapid rotation. The central density at the first core formation epoch is distributed for 
3 X 10^^ cm~^ ^ '^c ^ lO^^cm"^. In group A, models that fragment at relatively low den- 
sity have strong magnetic fields or rapid rotations, while models that fragment at relatively 
high density have weak magnetic fields and slow rotations. In addition, no fragmentation is 
expected for 10^^ cm^^ < nc < lO^^cm"^. This means that the core fragments just after the 
first core formation, and fragmentation never occurs if the core fails at this epoch. 

In group A, fragments have mutual separations of 3 — 300 AU. Since we stopped calcu- 
lation in group A after fragmentation occurs, we do not follow the cloud evolution until the 
protostar is formed {rtc — 10^^ cm~^). Thus, fragmentation may occur again for these models 
in the second collapse or protostellar phase. In that case, it is expected that hierarchical 
stellar systems are formed. 

Some models that do not show fragmentation in the adiabatic phase fragment in the 
second collapse phase (group B). The lower panel in Figure [TO] shows that the number of 
fragmentation models smoothly increases with the cloud density for > 10^^ cm~^. As 
shown in Figures [3] and [H the central region has a considerably we ak magnetic field in th e 



second collapse phase owing to the Ohmic dissipation (for detail, see lMachida et al. Il2007af ). 
Thus, the magnetic braking is not so efficient that the central region can spin up as the cloud 
collapses. Due to this accelerated rotation, fragmentation is induced in the second collapse 
phase. Fragments in group B have a typical separation of 0.01 — 0.3 AU. 

When the central density exceeds Uc > 10^^ cm~'^, the second core (or protostar) is 
formed. Even after the second core formation, fragmentation frequently occurs. In group C, 
fragmentation occurs just after the second core formation, similar to the fragmentation in 
the adiabatic phase. Fragments in group C have a typical separation of 0.005 — 0.07 AU. 

The fragmentation process both in the adiabatic and protostellar phases is thought 
to differ from that in the second collapse phase. Both in the adiabatic and protostellar 
phases, fragmentation occurs after the respective cores are formed. Fragmentation easily 
occurs in a quasi-static core, and the core has enough time to amplify the perturbation that 
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induces fragmentation. On the other hand, fragmentation occurs owing to cloud rotation 
in the second collapse phase. Fragmentation does not occur in the isothermal phase when 
the cloud collapses in a self-similar fashio n, because the self-similar solution with 7 = 1 
exists even in a rotating collapsing cloud (IMatsumoto et al. 1 119971 : iMatsumoto &: Hanawa 
19991 ). On the other hand, the clouds with 7 > 1 form rotating disks and stop the 



contraction, because n o self-similar s olutio n exists (ISaigo et al. 1 120001 ). These clouds can 
fragment as shown in ISaigo et al. I (120041 ). In addition, since the cloud collapses with 
7 = 1.1, the non-axisymmetr ic pe rturbation manages to evolve in this phase as shown in 



Hanawa fc Matsumoto I (I2OOOI ) and iLail (I2OOOI ). Thus, the fragmentation process in groups A 



and C is different from group B, which causes a gap between group A and group B. However 
we cannot distinguish group B from group C, because the number of fragmentation models 
smoothly increases after 11^ > 10^^ cm~^. As a result, the densities when fragments appear 
are divided into two groups: 10^^ cm~^ ^ '^c ^ 10^^ cm~^ and ric ^ 10^^ cm~'^. Histogram of 
separation shows that the separations between fragments are clearly divided into two groups 
(Fig. [To] right panel): models anticipating wide (3 — 300 AU) separations and those antici- 
pating narrow (0.005 — 0.3 AU) separations. However, this feature of two distinct groups of 
separation may be smoothed out as long as the protostar evolves. In addition, we stopped 
the calculation after the first fragmentation occurs or the protostar is formed. Thus, further 
fragmentation (or protostellar formation) may occur around a protostellar disk and the first 
core. Thus, although further lengthy calculations are needed to know the real distribution of 
binary separations, we may observe two distinct groups of binary in young clusters composed 
of extremely young stars. 



5.3. Does Magnetic Field Suppress or Promote Fragmentation? 

In Figure [TTl the fragmentation epochs are plotted against the ratio of the initial mag- 
netic energy to the rotational energy (_Emag/-Erot)- The right axis in this figure means the 
Jeans length that is calculated from the number density (right axis), in which the temperature 



is rela ted to the number density by the spherical symmetric calculation (IMasunaga fc Inutsuka 



2OOOI ). Only fragmentation models in the magnetized clouds are plotted in this figure. The 
figure shows that the fragmentation epochs shift to high density as £'mag/-E'rot increases. This 
indicates that cloud rotation promotes fragmentation, while the magnetic field suppresses 
fragmentation. This is valid through all the phases of cloud evolution from the isothermal 
to the protostellar phases. When the magnetic energy is lower than the rotational energy 
Emag < Ej-ot in the initial cloud, many models result in fragmentation in the adiabatic phase. 
On the other hand, when the magnetic energy is predominant over the rotational energy 
-^mag > Ej-ot, fragmentation occurs mainly in the second collapse and protostellar phases. 
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and some models do not show fragmentation. In this figure, there are some exceptions. Two 
models show fragmentation in the second collapse phase even if _Emag/-Erot < 0.1. In these 
models, fragmentation does not occur in the adiabatic phase, because the non-axisymmetric 
perturbation is moderately grown and an induced spiral pattern effectively removes the an- 
gular momentum from the center of the cloud. Thus, fragmentation is postponed. 

Observations have shown that molecular cloud c ores have high magnetic energies, but 
low rotational energies JCrutcheijEoggl : ICasellill2002h . which indicates ^'magZ-^rot ^ 1 Fig- 
ure [TT] indicates that these clouds barely fra gment. This conflicts with o bservations that the 



majority of stars are born as binary (e.g., iDuquennoy &: Mayor Ill99ll ). In this study, we 



calculated the evolution of a magnetized cloud in which the initial magnetic fleld lines are 
aligned to the rotation axis. When magneti c field lines are not align ed to the rotation axis, 
cloud evolution may be changed. However, iMachida et al. I (l2006a( ) showed that when the 
magnetic field lines are not aligned to the rotation axis at the initial state, the magnetic 
braking is more effective, and the angular momentum is effectively tr ansferred outwardly . 
This implies that fragmentation is more suppressed. On the other hand, iPrice fc Batd (120071 ) 
showed, in their ideal MHD calculations, that when initial cloud has a considerably distorted 
structure, fragmentation can occur even in a strongly magnetized cloud. 

To understand observational binary frequency, we may need to calculate the cloud evo- 
lution in a non-ideal regime, parameterizing the cloud shape and the angle between the 
rotation axis and the magnetic field lines, besides the non-axisymmetric perturbation, and 
the magnetic, rotational, and thermal energies. 



5.4. Fragmentation Condition and Wide/Close Binary Formation 

We investigated the cloud evolutions controlled by four parameters. Each parameter 
corresponds to the magnetic b, rotational u, thermal energies Uq, and the amphtude of the 
non-axisymmetric perturbation of the initial cloud. As shown in the previous sections, 
fragmentation significantly depends on the magnetic field strength and rotation rate of the 
initial cloud, while it slightly depends on the thermal energy and the amplitude of the non- 
axisymmetric perturbation. Thus, in Figure [121 we have plotted the final state of each cloud 
against the magnetic field strength and rotation rate of the initial cloud. In this figure, x- 
and y-axes are parameters b and u, respectively. Upper and right axes mean, respectively, 
the ratio of the magnetic (70) and rotational (/5o) energies to the gravitational energy of the 
initial cloud. Bottom axis means the mass to magnetic fiux ratio (M/<l>)norm normalized by 
the critical value (see, Eqs 10-12). In this figure, some models are located in the region of 
7o > 1 or (M/$)norm < I5 which indicates that the cloud is magnetically subcritical. In these 
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models, the initial cloud is magnetically subcritical as a whole, while the central region is 
magnetically supercrit ical since we adopted the Bonnor-Ebert density profile (for detail, see 



Machida et al. 



2007al ). Thus, any cloud can collapse in a sufficiently shorter timescale than 



the ambipolar diffusion timescale. 

In Figure [T2| three types of symbols (circle o, diamond o, and cross +) are plotted against 
the initial magnetic fields and angular velocities. Each symbol includes three models with dif- 
ferent {A^ = 0.01, 0.2, and 0.4). Symbols indicate that (i) Circle o: At least one model out 
of three shows fragmentation in the adiabatic phase (10^^ cm~^ ^ '^c ^ lO^^cm"^), (ii) Dia- 
mond o: in models that do not show fragmentation in the adiabatic phase [ric < 10^^ cm~^), 
at least one model out of three shows fragmentation in the second collapse or protostel- 
lar phase (n^ > 10^^ cm"'^), and (iii) Cross +: models show no fragmentation through any 
phases of cloud evolution (lO^cm"'^ < < 10^^ cm^'^). As shown in Figure [T^ we can 
divide the parameter space into three regions: wide fragmentation, close fragmentation, and 
non-fragmentation, which correspond to the formation of wide binary, close binary, and sin- 
gle star. The separations between fragments are clearly divided into two classes (wide and 
narrow separation) according to the fragmentation epoch, as shown in Figure [10 

Figure [12] indicates that models with strong magnetic fields and rapid rotations are 
distributed in the wide binary region. Models in the wide binary region can fragment in the 
adiabatic phase, and fragments have separations of 3 — 300 AU. However, not all models in 
this region fragment in the adiabatic phase. Even among models in this region, some models 
show fragmentation only in the second collapse or protostellar phase, and some models show 
no fragmentation through any phases of cloud evolution. For example, models {b, u) = 
(0.01, 0.05), which are located in the wide binary region in Figure [121 show fragmentation 
in the adiabatic phase for A^ = 0.01 (Fig. [1]) and 0.4 (Fig. [6]), while fragmentation occurs 
only in the second collapse phase for the model with = 0.2 (Fig. H]). In models (6, u) 
= (0.001, 0.07), fragmentation occurs in the adiabatic phase for the models with A^ = 0.01 
(Fig. [1]) and 0.2 (Fig. Hj), while the protostar is directly formed without fragmentation for 
the model with A^ = 0.4 (Fig. [6]). Thus, a small difference in the amplitude of the non- 
axisymmetric perturbation at the initial state has a possibility to induce a large difference 
in the fragmentation epoch. In some models located in the wide binary region, we stopped 
calculation before the protostar was formed when the fragmentation occurred in the adiabatic 
phase, because the fragments are likely to escape from the finest grid. Thus, fragmentation 
can occur again in each fragment in the second collapse or protostellar phase, which seems 
to lead to a hierarchical stellar system. 



Models distributed in the close binary region have a weaker magnetic field and slower 
rotation than models in the wide binary region. However, the tendency of rotation to promote 
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fragmentation, and the magnetic field to suppress fragmentation is the same. Models in the 
close binary region can fragment in the second collapse or protostellar phase, and fragments 
have separations of 0.005-0.3 AU as shown in Figure [TUl However, similarly to wide binary 
models, even in models in the close binary region, some models do not fragment in either 
the second collapse or protostellar phases. For example, in models (6, u) = (1, 0.05), 
fragmentation occurs in the protostellar phase for the model with = 0.01 (Fig. [1]), while 
protostars are directly formed without fragmentation for the models with = 0.1 (Fig. H]), 
and 0.4 (Fig.E]). 

Models in the single-star region have stronger magnetic fields and slower rotations than 
those in the wide and close binary regions. In all models in this region, the protostar is 
directly formed without fragmentation. As a result, a single compact sphere is formed at the 
center of the cloud, as shown in models {b, uj) = (4, 0.007) located in the lowest right corner 
of Figures [1], HI and [61 These compact spheres are expected to evolve into single stars. 

Models located in the wide/close binary regions have the possibility of forming the 
wide/close binary. As for fragmentation, it is difficult to forecast the cloud evolution from 
the initial conditions, and a small difference in the initial state sometimes results in different 
final outcomes. For example, even models located in the wide binary region have three 
possibilities: formation of a wide binary, a close binary, or a single star. That is, the 
parameters in wide/close binary regions are the necessary conditions, not sufficient conditions 
for wide/close binary formation. 



5.5. Comparison with Observation 



As shown in Figures [TT] and [121 our results indicate that cloud rotation promotes frag- 
mentation, but the magnetic field suppresses fragmentation. In this section, to quantify the 
magnetic field and cloud rotation in the molecular cloud core, we compare our results with 
observations. 

Both magnetic field and rotation rate are ob served in molecul ar cloud TMC-I C. TMC- 
IC has /5 = 1.2 X 10-3 JCoodman et al. lllQQsh . and 7 < 0.08 JCrutcheil Egooi ). As for 
magnetic field strength, only maximum value was estimated for observational limitation. 
When a cloud has f3 = 1.2 x 10"^, Figure [12] indicates that (i) fragmentation occurs in the 
adiabatic phase and the wide binary can form when 7 < 10"^, (ii) fragmentation occurs 
either in the second collapse or protostellar phases, and the close binary can form when 
10^'^ < 7 < 0.45, and (iii) only a single star forms without fragmentation when 7 > 0.45. 
Since observation showed 7 < 0.08, it is possible to form a binary system in TMC-IC. 
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In the observations of molecular cloud core, ICrutcherl (119991) showed that the mag netic 
energy i s con aparable to gravitational energy (7 ~ 0.5), while iGoodman et al. I (119931 ) and 



Casellil (|2002| ) showed that rotational is much smaller than gravitational energy {j3 ~ 0.02). 



Thus, these observations indicate that it is difficult to form a (wide) binary system as shown 
in §5.3. Except for TMC-IC, however, there are few clouds in which both magnetic field and 
rotation rate are observed. In addition, the observed magnetic and rotation energies might 
not reflect the distribution of these energies in the majority of cloud cores, because we can 
only observe a limited range of the magnetic field strength and rotation speed: we cannot 
observe weak magnetic field strengths and slow rotation speeds. Furthermore, we have to 
measure the magnetic field strength and rotation rate at the same radius or the same number 
density in the molecular cloud core. To determine magnetic field strength and rotation rate 
precisely, we need a high-resolution observational facility such as the future ALMA. 



5.6. Comparison with Previous Work 



Fragmenta tion conditions i n the adiab atic phase (i.e . , the wide binary region) were 
investigated by iMachida et al. I (l2005bl ). and lPrice &: Bate! (120071 ) using the ideal MHD ap- 



proximation, which is not valid in the high density gas region with 10^^ ^ '^c ^ 10^^ cm~^. 
In this paper, we studied cloud evolution in the non-ideal MHD regime. However , the 
wide binary region in Figure [T2] corresponds well to results of iMachida et al. I (l2005bl ) . As 
shown in Figure [TOl many models located in the wide binary region fragment in the range 
of 10^^ cm~^ ^ '^c ^ 10^^ cm~'^. In addition, since the density outside a small central part of 
the cloud is < 10^^ cm~^, the Ohmic dissipation is not so effective at the fragmentation 
epoch. From this, there must be only a little difference between the ideal MHD and non-ideal 
MHD calculations for fragmentation appearing in the adiabatic phase. 



Price &: Batd (120071 ) investigated the evolution of magnetized clouds and fragmentation 
conditions with their MHD-SPH code. Their results are qualitatively the same as ours: 
the magnetic field suppresses fragmentation. However, there is a small difference. They 
adopted a rotatio n energy of /3n = 0.0 05 for the initial clouds. Fragmentation occurs when 
(M/$)norm > 5 iu lPrlce fc Batd (120071 ). while fragmentation occurs when (M/ $)norm > 10 in 
our ca lculations for models with /3o = 0.005, as shown in Figure [121 Thus, in iPrice fc Bate 



( 120071 ). fragmentation occurs in stronger magnetic fields than ours. 



Fragmentation is promoted by rotation, but it is suppressed by the magnetic field. 
In strongly magnetized clouds, since the cloud rotation is removed from the central re- 
gion by the magnetic braking and outflow, fragmentation rarely occurs. Outflow is closely 
related to the magnetic braking: both outflow and magnetic braking are caused by the tor- 
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sional Alfven wave or the magnetic tension force generated by rotation of the central core. 
Thus, the outflow is one of the proofs that the Alfven wave is properly resolved. When 
the Alfven waves are resolved, the outflow always appears after the adiabatic core ( or first 
core) is formed. In the collapsing clouds, no outflow appears in iPrice fc Bate while 
the outflow appears in t he gr id -based MHD simulations, su c h as Tomisaka J 19981. 2002 ) , 



Matsumoto k Tomisakal (12004 ). iMachida et al. I (12004 l2005bf ) . iFromang et al. I (120061 1. and 



Banerjee fc Pudrita (120061 ). It is expected that the angular momentum transfer caused by 



the torsional Alfven wave (or the magnetic tension force) may not be correctly resolved in 



Price fc Batd (120071 ). thus excess angular momentum around the center of the cloud is left. 



Therefore, in their calculation, fragmentation occurs more frequently than ours. 



6. SUMMARY 

We investigate the cloud evolutions for 10^ cm"'^ < ""-c ^ 10^^ cm~^ in a large parameter 
space. In this study, we systematically calculated 147 models with different magnetic field, 
rotation, thermal energies, and the amplitude of the non-axisymmetric perturbation of the 
initial cloud. These calculations indicate that 

1. Fragmentation significantly depends on the magnetic field and rotation, but slightly de- 
pends on the thermal energy and the amplitude of the non-axisymmetric perturbation 
of the initial cloud. 

2. The magnetic field suppresses fragmentation, and rotation promotes fragmentation 
through all phases of cloud evolution. 

3. The distribution of the separations between fragments are clearly divided into two 
classes: fragments formed in the adiabatic phase have wide separations as 3 — 300 AU, 
while fragments formed in the second collapse and protostellar phases have narrow 
separations as 0.007 — 0.3 AU. 
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Table 1: Model parameters of the magnetic fields 
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Table 2: Model parameters of the angular velocities 
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Fig. 1. — Final states for models with = 0.01 against parameters b, and u. Density 
(color-scale) on the cross-section of the z=0 plane are plotted in each panel. The grid level, 
maximum number density (nmax), and grid scale are denoted inside each panel. Background 
colors indicate that Blue: fragmentation occurs in the adiabatic phase, Red: fragmentation 
occurs only in the second collapse or protostellar phase. Gray: no fragmentation occurs 
through all phases of the cloud evolution, and Green: the cloud no longer collapses. 
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Fig. 2. — Cloud evolution for model (6, cj, A^) = (0.1, 0.05, 0.01). Density (color-scale) and 
velocity distribution (arrows) on the cross-section of the z = plane (upper panels), y = 
plane (middle panels), and x = plane (lower panels) are plotted. Panels (a) through (f) are 
snapshots at the stages (a) ric = 1.7 x lO'^cm"^ (/ = 3 — 5; initial state), (b) 7.3 x 10^^ cm"'^ 
(/ = 16-19; adiabatic phase), (c) 3.6 x 10^'' cm^'^ (/ = 18-20; second collapse phase), (d) 
3.5 X 10^2 cm-3 (^i ^ 25, 26; protostellar phase), (e) 7.0 x 10^^ ^^-3 ^ 26, 27; protostellar 
phase), and (f) 5.2 x 10^^ cm~^ (/ = 26, 27; calculation end), where I denotes the level 
of the subgrid. The red lines in panels b and c indicate the first core. The black lines in 
panels d-/ indicate the protostar (or second core). The thick orange lines in panel b indicate 
the border between the infalling and outflowing gases (contour of = 0). The level of the 
subgrid is shown in the upper left corner of each upper panel. The elapsed time t, maximum 
number density nmax, and arrow scale are denoted in each panel. 
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Fig. 3. — Evolution of (a) oblateness, (b) axis ratio, (c) angular velocity normalized by the 
square root of the central density fic/(47rG'pc)^^^) and (d) magnetic field strength normalized 
by the square root of the central density B^c/ (47rCg gPc)^^^- All quantities are plotted against 
the central density ric for model {b, ur, A^) = (0.1, 0.05, 0.01). 
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Fig. 5. — Cloud structure for model {b, u, A^) = (0.01, 0.05, 0.2) at the end of the 
calculation with different scales. Density (color-scale) and velocity distribution (arrows) on 
the cross-section of the z = plane (upper panels), y = plane (middle panels), and x = 
plane (lower panels) are plotted. The red dotted line in panel b indicates the first core. The 
thick orange lines in panel b indicate the borders between the infalling and outflowing gases 
(contour of Vz = 0). The level of the subgrid is shown in the upper left corner of each upper 
panel. Arrow scale is also denoted in each panel. 




Fig. 7. — Same as Fig. [5l Final state for model {b, u, A^) = (0.01, 0.07, 0.4) is plotted 
with different grid scales. 



-41 - 




Fig. 8. — Evolution of models with different ao- Panel (a): Oblatencss eob, (b): axis ratio 
^ar, (c) Magnetic field normalized by the square root of the central density B^c/ (47rCg qPc)^^^, 
and (d) angular velocity normalized by the square root of the central density flc/ (47rGpc)^^^) 
all of which are plotted against the central density ric- 
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Fig. 9. — Cloud structures around the center of the cloud for models with different ao 
(ao = 0.1,0.3,0.5, and 0.7) at the fragmentation epoch. Density (color-scale) and velocity 
distribution (arrows) on the cross-section in the z = plane are plotted. The elapsed time 
t, maximum number density rimax, and arrow scale are listed in each panel. The level of the 
subgrid is shown in the upper left corner of each panel. 
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Fig. 10. — Central panel: Separations between furthermost fragments as a function of 
the density at the end of the calculation. The Jeans length evaluated from the spherical 
symmetric calculations is also illustrated by the dotted line. The results of all models where 
fragmentation is observed are plotted. Models are classified into three groups according to the 
fragmentation epoch: fragmentation in the adiabatic phase (A), second collapse phase (B), 
and protostellar phase (C). Right panel: Histogram of fragmentation models as a function 
of the separation r. Bottom panel: Histogram of fragmentation models as a function of the 
fragmentation epoch measured by the central density Uc- 
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Fig. 11. — Number density at the fragmentation epoch as a function of the ratio of the 
magnetic to the rotation energies £'mag/-E'rot of the initial cloud. Right axis indicates the 
Jeans length. The dashed hues mean the borders between different evolutional phases. 
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Fig. 12. — Fragmentation condition against parameters b and u. Symbols mean that circle o: 
fragmentation occurs in the adiabatic phase, diamond o: fragmentation occurs in the second 
collapse phase, cross +: no fragmentation occurs through all phases of the cloud evolution. 
Models in the wide and close binary regions have separations in the range of 3 — 300 AU 
(wide binary region ), and 0.007 — 0.3 AU (close binary region). Models in the single-star 
region do not show fragmentation in any evolution phase. Upper and right axis indicate the 
ratio of the magnetic 7 (upper) and rotational (3 (right) energies to the gravitational energy. 
Bottom axis means the mass-to- magnetic flux ratio M/$ normalized by the critical value 
(see, Eq[l2]). 



